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Oscillating fields can make domain patterns change into various types of structures. Numerical 
■ simulations show that concentric-ring domain patterns centered at a strong defect are observed under 

a rapidly oscillating field in some cases. The concentric-ring pattern appears near the threshold of 
spatially-uniform patterns in high-frequency cases. The threshold is theoretically estimated and 
the theoretical threshold is in good agreement with numerical one in a high-frequency region. The 
theoretical analysis gives also good estimations of several characteristics of domain patterns for 
high-frequencies. 
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PACS numbers: 89.75.Kd, 75.70.Kw, 47.20.Lz, 47.54.- 
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Rapidly oscillating fields cause interesting phenomena in a wide variety of systems. Those phenomena are often 
discussed in the view of stabilization of an unstable state. One of the simplest examples is the problem of Kapitza's 
inverted pendulum, which was generalized by Landau and Lifshitz The key idea is to separate the dynamics 
^ ' into a rapidly oscillating part and a slowly varying part. The method has been applied to various systems, e.g., the 
classical and quantum dynamics in periodically driven systems 0, the stabilization of a matter-wave soliton in 
two-dimensional Bose- Einstein condensates without an external trap [J, [H, Q , and magnetic domain patterns traveling 
at a slow velocity under a rapidly oscillating field [7[. 
, Domain patterns are observed in a wide variety of systems, and they show many kinds of structures (see, for 
example, Refs. p| H| [l(| and references therein). Magnetic domain patterns are one of their good examples. While 
they usually exhibit a labyrinth structure under zero field, they also show other kinds of structures under an oscillating 
field. For instance, parallel-stripe and several kinds of lattice structures are observed in experiments and numerical 
simulations [TTI . 12|, EH . Moreover, traveling patterns 0] and more interesting patterns, e.g., spirals and concentric-ring 



CNJ ■ patterns [1J, [15| , have been observed, depending on the strength and frequency of the field. Spirals and concentric- 
ring patterns appear under a large-amplitude and high-frequency field in experiments, and the field range where they 
. appear is not wide. 

In this paper, we investigate effects of an oscillating field by numerical simulations and theoretical analysis, focusing 
on the emergence of concentric-ring magnetic domain patterns surrounding a strong defect. Recently, two theoretical 
methods were proposed to investigate the effects of an oscillating field on pattern formation in ferromagnetic thin 
films Q ■ One gives the "time-averaged model" and the other gives the "phase-shifted model" : The former is derived 
by averaging out rapidly oscillating terms, and the latter includes the delay of the response to the oscillating field. In 
this paper, the time-averaged model is applied to theoretical analysis, since it is suitable for discussing "stationary" 
domain patterns which oscillate periodically but are unchanged in terms of a long-time average. The theoretical 
line of the threshold for nonuniform patterns is derived from the time-averaged model. The theoretical threshold is 
consistent with the numerically simulated one in a high-frequency region, where a concentric-ring pattern appears 
;J] ■ around the defect. 

In fact, there are several techniques to study domain patterns under a rapidly oscillating field theoretically. Applying 
a multi-time-scale technique [HI Ht| . one can obtain more complex equations in a better approximation than the 
time-averaged model. In other words, the time-averaged model corresponds to the lowest orders of multi-time-scale 
expansions. In this paper, the time-averaged model is employed since it has a simple form and is efficient enough to 
discuss the appearance of concentric-ring patterns in a high-frequency region. 

The creation of a concentric-ring pattern can have different mechanisms. One of them is boundary conditions, and 
the strong defect is a kind of boundary condition. The selection of a pattern depends on boundary conditions as well as 
the field frequency or other parameters [l8| . For example, in nematic liquid crystals under a rotating magnetic field, it 
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is sometimes observed that the center of concentric rings nucleated by a dust particle moves away from it [19( ; Faraday 
experiments of viscous fluid and granular layers in round cells show concentric- ring patterns or spiral patterns [2(J 
in some cases. On the other hand, concentric-ring patterns can also appear spontaneously. In fact, spiral patterns 
as well as concentric-ring patterns appear in the absence of a strong defect under some conditions [1J,|22|. However, 
we will not consider spontaneously created concentric-ring patterns in this paper since those patterns are beyond the 
scope of this paper. 

The rest of this paper is organized as follows. In Sec. [HI the model of our system is introduced and numerical results, 
i.e., the phase diagram for concentric-ring patterns and profiles of the domain patterns, are exhibited. In Sec. IIII[ we 
discuss the threshold for nonuniform patterns, employing the time-averaged model. Moreover, several characteristics 
of a domain pattern estimated from the time-averaged model are compared with those from the numerical simulations. 
The mechanism for the appearance of a concentric-ring pattern is discussed in Sec. lIVI Conclusions are given in Sec.lVl 



II. NUMERICAL SIMULATION 

Our model is a simple two-dimensional model (see Refs. @, HH [HI j an d references therein). The Hamiltonian of the 
model consists of four energy terms: Uniaxial anisotropy energy H an [, exchange interactions Hj, dipolar interactions 
-ffdb and the interactions with the external field H CK . We consider a scalar field 4>(r), where r — (x, y). The anisotropy 
energy is given by 

^ (i) 

where A(r) is employed to express the effect of defects or the roughness of a sample. This term implies that the values 
4>{r) = ±1 are preferable. The positive and negative values of <j>{r) correspond to up and down spins, respectively. 
The exchange and dipolar interactions are described by 

Hj = sj ir mpi (2) 

and 

i?di = 7 J drdr'(p(r)cp(r')G(r,r'), (3) 

respectively. Here, G(r,r') ~ \r — r'\~ 3 at long distances. These two terms are competing interactions: Hj implies 
that <j){r) tends to have the same value as neighbors, while iJ<ji implies that <p(r) prefers to have the opposite sign to 
ones at some distances. The interactions with the external field is given by 

H cx = -h(t) J dr<f>(r). (4) 

Here, we consider a spatially homogeneous and rapidly oscillating field, 

h{t) = hosinwt. (5) 

From Eqs. ([T])-([l]), the dynamical equation of the model is described by 

d<t>{r) = SjH ani + Hj + H di + H CK ) 
dt 6cf>{r) 

= oA(r)[^(r) - <Mr) 3 ] + /3V 2 0(r) - 7 J dr'0(r')G(r, r') + hit). (6) 

The numerical procedures are almost the same as those of Refs. @, [H, [24|. For time evolution, a semi-implicit 
method is employed: The exact solutions and the second-order Runge-Kutta method are used for the linear and 
nonlinear terms, respectively. For a better spatial resolution, a pseudo-spectral method is applied. Namely, the time 
evolutions are calculated for the equation in Fourier space corresponding to Eq. ((6|) , 

^ = aU - f)\] k - i(3k 2 + jG k )cb k + h(t)S h , (7) 
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FIG. 1: (Color online) Phase diagram for the concentric-ring patterns surrounding a strong defect at the center. The horizontal 
and vertical axes are the frequency (uj) and the amplitude (ho) of the external field, respectively. The red solid lines and blue 
dashed line are obtained from the numerical simulations and theoretical analysis, respectively. The values of u> and ho of each 
snapshot are as follows: (a) oj/2-k = 2/3 ~ 0.667 and ho = 3.05, (b) uj/2h = 2/3 and ho = 2.8, (c) lu/2tt = 0.1 and ho = 1.2, 
(d) u)/2n — 0.5 and ho = 2.5, and (e) uj/2h = 0.8 and ho = 3.56. 



where [•] denotes the convolution sum and Gk is the Fourier transform of G(r,0). Here, we define G(r,0) = l/|r| 3 . 
Then, one has 

Gk = aa - aik, (8) 

where k = \k\ and 

f°° dr 

a Q = 2Tr — , at = 2tt. (9) 
Jd r 2 

Here, d is the cutoff length of the dipolar interactions. In the simulations, we set d — ir/2, which results in ao = 4. 

The effect of a strong defect is incorporated in the anisotropy term, Eq. (Q}. Here, we put the strong defect at the 
center, i.e., the origin (r = 0), as follows: 

This condition implies that the spin at the center will not flip unless the applied field is too strong. The parameters 
in Eqs. (HJ)-© are given as a = 2, j3 — 2, and 7 = 2(3 /ai = 2/ir. The simulations are performed on a 128 x 128 lattice 
with periodic boundary conditions. 

The concentric-ring patterns simulated by the numerical calculations appear in a limited region of the frequency (u) 
and the amplitude (ho) of the external field. Figure [1] shows the u>-ho phase diagram for concentric-ring patterns. The 
solid lines and the dashed line are drawn by using the results of the numerical simulations and theoretical analysis, 
respectively. The theoretical analysis is explained in Sec. IIIIl Concentric-ring patterns are seen only in the region 
between the upper and lower solid lines. Above the upper solid line, one sees only spatially-homogeneous patterns 
except for the vicinity of the center. On the contrary, stripes, labyrinth or lattice structures appear below the lower 
solid line. 

Actually, it is difficult to find the exact boundaries of the region where concentric-ring patterns appear. One 
may think that concentric-ring patterns can appear right on the numerical threshold line in a low-frequency region. 
However, even if a concentric-ring pattern could appear on the threshold, it would be hard to find its exact value. In 
a low-frequency region, the boundary between spatially- homogeneous patterns and nonuniform patterns (e.g., stripes, 
labyrinth, and lattice structures) is sharp. In contrast, in a high-frequency region, the boundaries between concentric- 
ring patterns and other patterns are unclear: They are crossover lines rather than transition ones. Actually, for high 
frequencies, a few concentric rings around the strong defect coexist with other patterns (i.e., stripes, labyrinth or 
lattice patterns) in the region near the boundary lines of the concentric-ring-pattern region. In other words, a few 
concentric rings start to appear at the lower boundary, and the number of rings grows as ho increases. At the upper 
boundary, nonuniform patterns disappear except for the vicinity of the strong defect. 
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FIG. 2: Profiles of the domain patterns corresponding to the snapshots (a) and (b) in Fig. \T\ (a) uj/2n — 2/3 ~ 0.667 and 
ho = 3.05, (b) oj/2-k — 2/3 and ho = 2.8. The upper and lower thin curves in each figure are the profiles at t = (n + 1/2)T and 
t = nT, respectively, where n is an integer and T = 2n/ui. The middle thick curve is the time-averaged profile. 

The profile of each domain pattern is useful to see the time dependence of the pattern. The profiles corresponding 
to the snapshots (a) and (b) in Fig. [T] are shown in Fig. [5] The profile is a section which includes the center (the 
defect) and is perpendicular to x axis (the horizontal axis). The profiles show that the pattern is oscillating without 
deformation except for the vicinity of the center. The amplitudes of the oscillation and that of the pattern (i.e., the 
difference between the maximum and minimum values of <j>(r) except for the vicinity of the defect at a certain time) 
depend on the amplitude and frequency of the field. 



III. THEORETICAL ANALYSIS 

Since a rapidly oscillating field is applied, Kapitza's idea [l[ is applicable to the analysis of the pattern formation 
in this system. In fact, the profiles in Fig. [2] validate the use of the idea. In other words, the results in Fig. [5] justify 
the fact that the variable <f>(r,t) in Eq. ([6]) consists of a spatially homogeneous oscillating term <fio(t) and a slowly 
varying term $(r,i). In this section, the domain patterns in the absence of a defect are discussed by employing the 
time-averaged model @- Namely, A(r) in Eq. (fT0|) is replaced by unity, i.e., A(r) = 1. 

First of all, let us consider the spatially homogeneous oscillating solution 4>o{t) of Eq. ((6]). Then, one has 

0o - a(0 o - 0o) ~ 70o J dr'|r'r 3 + h(t). (11) 
Its solution can be approximately written as 

0o =psin(wi + 8), (12) 

where 8 is a phase shift which comes from the delay of the response to the field. Substituting Eq. (p~2|) into Eq. (fTTj) 
and omitting high-order harmonics (i.e., sin3a->£), one can evaluate p and 8. The value of p is obtained from the 
following equations [7|: 

^aV-^70P 4 + (- 2 +%V = ^, (13) 

where i]q = a — 70,0- 
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Now, we consider the equation for the slowly varying variable $(r, t) which describes a spatially dependent pattern. 
Substituting </>(r, t) — <fio(t) + $(r, t) into Eq. ^ and averaging out the rapid oscillation, one obtains the time-averaged 
model 



d<f>(r] 
dt 



= (l - ^/A a$(r) + /3V 2 $(r) - 1 J dr' ^'J - a$ 



(rf 



(14) 



The linear stability of Eq. (TTJ} leads to the theoretical curve in Fig. Q] corresponding to the threshold for the 
existence of nonuniform patterns. Substituting $(r) = J^k exp(ik ■ r)$fc into Eq. (Tl"4"|) , one has the linear part of the 
equation written as 



dt 



(15) 



Here. 



Vk = ( 1 - )a — /3{k-ko) + f3k - ^a Q , 



(16) 



where k = \k\ and fc = aij/(2/3). Since % has the maximum value at k = k , the value of p for = gives the 

instability threshold p ci 



3a 



(a + (3kl - ja ) 



1/2 



(17) 



When p > p c , f]k is negative for all values of fe. In other words, the homogeneous pattern, i.e., 5>(r) = 0, is stable and 
no inhomogeneous pattern tends to appear for p > p c . The threshold curve for nonuniform patterns in Fig. [T] (the 
dashed line) is given by Eq. lfl3|) with p = p c . 

Now, let us discuss how a nonuniform pattern disappears near the threshold. Taking <I>(r,i) = Acosfco^ which is 
one of the simplest stable patterns and substituting it into Eq. (|14[) , we have 



ml 



7O 



3 ,2 

—uA 

4 



A cos kox — —A^ cos = 0. 



Neglecting the higher harmonics (i.e., cosSfco^); we obtain 
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(18) 



(19) 



The amplitude A of the pattern decreases monotonically in terms of p and vanishes at p = p c . Namely, the amplitude 
of the pattern diminishes near the threshold. This behavior is found in Fig. [21 The same behavior can be derived for 
a concentric-ring pattern, which needs more complex calculations. 

The validity of the above discussion is examined by comparing numerical results with theoretical estimates. Actually, 
Fig. [1] indicates that the theoretical threshold is in good agreement with the numerical one for high frequencies 
(oj/2tt > 0.5). More quantitative comparisons are given in Table Q] for w/27r > 0.5. The numerical and theoretical 
values of p and A are compared in it for the data near the threshold. They are obtained from the profiles of domain 
patterns. Namely, the one-cycle time sequence of the profiles is used in order to estimate them. The theoretical value 
of p is calculated from Eq. (fT3")) , and that of A is calculated in two ways: The value of p in Eq. (fT^|) is given by (I) the 
value from simulations and (II) the theoretical value. The values of A estimated from simulations and Theory (I) are 
in good agreement. Moreover, when the numerical and theoretical values of p are close, the value of A from Theory 
(II) also has a similar value to the corresponding numerical A. 



IV. DISCUSSION 



Now, we consider how concentric-ring patterns appear around a strong defect in this system, employing the time- 
averaged model. Actually, the time-averaged model (fi"4j) without a strong defect can also produce concentric-ring 
patterns. If the initial condition is given as <i>(r = 0) = 1 at the center and <i>(r) = —1 anywhere else, then 
concentric-ring patterns, as shown in Fig. [31 are demonstrated by the time-averaged model without a strong defect. 
The values of p for Figs. 03a) and EHb) are given by p — 0.69 and p — 5, which correspond to those for Figs. [Ha) 
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TABLE I: Numerical and theoretical values of p and A. The numerical value of p is obtained from the one-cycle sequence of 
the profiles of domain patterns. The theoretical value of p is calculated from Eq. (|13[1 . and the values of A for Theory (I) and 
Theory (II) are given by Eq. (|19p with the numerical and theoretical values of p, respectively. 




FIG. 3: Domain patterns simulated by using the time-averaged model (|14|l without a strong defect: (a) p — 0.69, and p = 0.5. 
The parameters for (a) and (b) correspond to those of Figs. [TJ a) and[TJb), respectively. 

andQJb), respectively. One sees that Figs.[3fa) andQJa) look very similar. They are observed near the threshold for 
nonuniform patterns, where the time-averaged model is valid. Near the threshold, the linear growth rate r\k given by 
Eq. (|16|) is very small, which means that the time scale of pattern formation is very slow. In fact, it takes more than 
10 times longer time to obtain Fig. [3Ja) than Fig. [3Jb). If the initial pattern were completely uniform, nonuniform 
patterns could not appear easily. In other words, the initial condition employed here causes concentric rings on a 
uniform initial pattern (except for the center) as one of the simplest patterns. This situation is very similar to that 
of Fig. [Ha): The strong defect in the original model behaves as a kind of boundary condition in the almost uniform 
pattern near the threshold. Since the growth rate of the instability is very slow near the threshold, no patterns except 
for concentric rings can grow for Figs. [3ja) orQJa). The growth rate is expected to be relevant to the correlation 
length, which is related to the number of concentric rings. From analogy with critical phenomena, one can estimate 

— 1/2 

that the correlation length would be proportional to r] ko , where i]k = r]k=k ■ Therefore, the correlation length 
diverges near the threshold, as rjk — > 0. This is the mechanism of the emergence of a concentric-ring pattern around 
a strong defect and also the reason why the boundaries of the concentric-ring-pattern region are unclear in the phase 
diagram in Fig. [1] 

In contrast, Fig.[3{b) has a few concentric rings, although stripes (or mazes) spread outside the rings. The stripes 
grow independently from the concentric rings because of rather large 77^ 's. The difference between Figs.[3jb) and[lfb) 
comes out not only because the time-averaged model is invalid far from the threshold, but also because the initial 
condition employed here is pretty different from the situation for Fig. [T^b). 

The time-averaged model is also invalid for low frequencies. This is evident from Fig. [T] in which the theoretical 
threshold line is not consistent with the numerical one for uj/2it < 0.5. The failure of the time-averaged model in a 
low- frequency region comes from the assumption </>(r, t) = <fio + $(r, t). The assumption is valid when the time scales 
of the rapidly oscillating part and the slowly varying part are well separated. For low frequencies, their separation is 
not sufficient. If a multi-time-scale technique were applied, a better theoretical threshold line would be obtained. 

The concentric-ring patterns shown in this paper are purely numerical results, although they suggest a possible 
mechanism about the formation of the patterns. Since the characteristics of domain patterns strongly depends on 
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experimental conditions and samples, it is rather hard to compare experimental data with the results obtained in 
this paper. However, we can suggest that concentric-ring patterns appear in a certain range of the field strength and 
frequency, which is located above a labyrinth-pattern region. In fact, also in experiments, spirals and concentric rings 
are often observed near the threshold of nonuniform patterns [13, EH, although they are not always centered at a 
strong defect but often move around. Incidentally, for typical ferromagnetic garnet films, the order of the characteristic 
domain width is about 10-100 /jm and that of the frequency for lattices, spirals, or concentric rings is about 0.1-100 
kH 0[H. 



V. CONCLUSIONS 



In this paper, effects of an oscillating field have been investigated, and the emergence of concentric-ring patterns 
surrounding a strong defect has been discussed. The numerical simulations show that the concentric-ring pattern 
appears in the high-frequency region near the threshold for nonuniform patterns. The simulated profiles of the 
concentric-ring pattern indicate that the pattern consists of two parts (except for the vicinity of the defect) : A 
rapidly oscillating spatially homogeneous part and a nonuniform pattern part. This fact assures that the time- 
averaged model is suitable for the theoretical analysis. The theoretical threshold line is in good agreement with the 
numerical one in a high-frequency region. When the rapidly oscillating field makes the state close to the threshold, 
the concentric-ring pattern appears due to the strong defect which is an effective boundary condition. 

In conclusion, the validity of the time-averaged model has been demonstrated in the presence of a rapidly oscillating 
filed. It gives the good estimate of the threshold for nonuniform patterns when the field frequency is high. Moreover, 
it is revealed that ideal and interesting patterns such as concentric-ring patterns can appear near the threshold, 
depending on boundary conditions. 
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